On the local existence of maximal slicings in 
spherically symmetric spacetimes 



Isabel Cordero-Carrion^, Jose Maria Ibanez^ and Juan Antonio 
O I Morales-Lladosa^ 

^ Departamento de Astronomi'a y Astrofi'sica, Universidad de Valencia, C/ Dr. Moliner 50, 
E-46100 Burjassot, Valencia, Spain 

E-mail: isabel.cordero@uv.es, jose .m. ibeaiezSuv . es , eoitonio .morales@uv. es 



Abstract. In this talk we show that any spherically symmetric spacetime admits locally a 
maximal spacelike slicing. The above condition is reduced to solve a decoupled system of first 
O , order quasi-linear partial differential equations. The solution may be accomplished analytical 

or numerically. We provide a general procedure to construct such maximal slicings. 

^ ' 1. Introduction 

'nI" , A maximal hypersurface is one that has vanishing mean extrinsic curvature, K = 0, K being 

the trace of the extrinsic curvature of the hypersurface. The name comes from the fact that the 
induced volume functional reaches a local maximum with respect the variations that keep fixed a 
given boundary. Maximal hypersurfaces were considered by Lichnerowicz [I] to solve Einstein's 
(3 ■ constraint equations, giving motivation for subsequent studies on the subject (see, for example 

. El El HIE])- In fact, the existence of maximal hypersurfaces is extensively used in Mathematical 

Relativity. This property is a very simple geometric assumption to establish general results for 
broad classes of spacetimes, for instance, local or asymptotically stationary or conformally flat 
^ ', spacetimes. 

$H ' We will use the term maximal slicing when referring to a (non intersecting) family of spacelike 

maximal hypersurfaces which locally foliates a certain domain of spacetime. This type of slicing 
has very nice properties as, for example: i) the well-known singularity avoidance capability 
[6], ii) it is well adapted to the propagation of gravitational waves [71 [8], and, iii) it gives the 
natural Newtonian analogous when, in addition, conformal flatness is imposed on each slice 
[9]. Maximal slicing condition has been recently used in the Fully Constrained Formulation of 
Einstein equations derived by the Meudon group [lOl [11] . 

In spite of their extended use, the existence of maximal slicings in spherically symmetric 
spacetimes (SSSTs) has been only established for vacuum and for some particular energy contents 
(see [Hi [121 [121 [12 HSl [Ml HZj). There is, as far as we know, no theorem stating that always it 
is possible to build a maximal slicings in a SSST. 

In this work, we aim to prove the local existence of maximal slicings in any SSST. We 
will follow a purely geometrical approach, independent of Einstein equations, according to 
complementary to the standard time evolution strategy [9l [121 [13 [13 [161 [T7]. 

Although our study is independent of the field equations, one by-product of our approach, 
which could be of interest in the field of Numerical Relativity, is that it provides a means to 



assess complex and sophisticated 3D numerical codes built to solve Einstein equations. 



2. Local existence of maximal slicings 

In this section we establish the following result: 

Theorem. Any spherically symmetric spacetime can be locally sliced by a family of maximal 
spacelike hyper surf aces. 

In order to prove the previous theorem, we derive a decoupled system of three first order 
partial differential equations that proves the local existence of a maximal slicings, and provides 
a general procedure allowing its construction. 

Let us start with the canonical form of the metric of a SSST, 

ds^ = Adt^ +2C dt dr + Bdr'^ + Dd^'^, (1) 

where = dO"^ + sin^ Odtp^ is the metric of the 2-sphere, A,B,C,D are smooth functions of 
t and r, and AB — < to ensure the Lorentzian character of the metric. In addition we 
choose the signature (— , +, +, +), and accordingly L) > 0. Partial derivatives with respect to r 

df df 
will be denoted as — — = /', and with respect to t as — = /. The spatial metric induced on 

or dt 
the hypersurfaces S^, defined by t = constant, is 7jj = disig{B , D , D s\v? 9) , where B > since 

we are considering spacelike hypersurfaces. Let n be the future pointing timelike unit normal 

to the hypersurfaces St, 



''=a[d-t-Bd-r)^ « = V^-^- 

The mean extrinsic curvature K of is related with the expansion of n, K = —V • n, where V 
is the covariant derivative with respect to the spacetime metric. In the given metric in Eq. ([T]), 
this relation is 

In the following, we assume that A,B,C and D are known functions. We look for a change of 
coordinates {t = i{t, r),f = f(t, r),9, y?} such that the hypersurfaces t = constant are maximal. 
We introduce two fields, X and Y, satisfying the commutation relation [X, y] = 0. This 
condition assures the existence of two coordinate parameters, namely i = i{t,r) and r = r{t,r), 
such that 

X = l,Y = -l-. (4) 

dt dr ^ ' 

Then, we decompose these fields asY = \Y, X = oY + bY^, with ^ 1, y • y-^ = 0, 6 / 0, 
and A > 0. The condition y^ = 1 is equivalent to 

^ = ^1 + ^1^' P = B''{-fC + e./WT^)^ (5) 

being / an unknown function to be determined, e = ±1 and l"^ = —AB + > 0. Fixing the 
d 

coefficient of — in the decomposition of y-*-, Y ■ Y-^ = leads to 
dt 

y^ = i+Ql^Q = i-c + -y^iL=\ . (6) 

dt ^dr' ^ \ ^Bj 



Consequently, the resulting fields are 



X = (a/0)| + (<.P + 6Q)A, r = A(/| + p|).,(y„„+|v7^|), (7) 

and ([5]) . The condition [X, y] = is then equivalent 



where we have taken into account Eqs. 
to 



' p' 


1 


[/] 




'aP + bQ' 


1 


[a/ + 61 


bp 




bp 


1 


bp A 




bp A 



(8) 



where p = P — fQ 



Now, we denote with K the trace of the extrinsic curvature of the new hypersurfaces 
t = constant. The condition K = and the commutation relation provide 3 equations for 
4 unknown functions, a, b, /, A. Taking into account that L> is a scalar under the above change 
of coordinates, we can add, without loss of generality, the following coordinate condition 



(9) 



saying that the metric on the hypersurfaces t = constant is written in isotropic conformally flat 
form. From Eqs. ([3]) and ([9|), the condition K = is equivalent to 



2Y (X-Y)- 3X (y2) + 



4 Y (y2) 



X • y = 0. 



Prom the decompositions d?]), Eqs. @ and (|lUp are expressed as A 
viewed as a definition of r in terms of A, and as 



(10) 

D/r, that can be 
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After some algebraic calculations, the previous definition of A and Eqs. 
equivalent to 

.^/a + qv D + QD' \ 

-=b^[—x 
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and 



fUp'- = p'-p^l + f-fi 
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and (jlip are 
(12) 

(13) 
(14) 
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Notice that Eq. ()15p involves only / when p and Q are written explicitly in terms of /. 

First, Eq. (fT5]) can be solved for /. Second, Eqs. (fT3]l and (fH|l can be solved for A and 
b. Finally, a can be obtained from Eq. (|12p . Assuming that A,B,C,D are continuously 
differentiable functions, the initial value problem with respect to this set of equations has always 
local (both in space and time) solution [18] (which is also continuously differentiable). Therefore, 
we have proved the announced theorem. 



Notice that in order to solve this set of equations, it can be useful to distinguish two different 
cases, / = and f ^ 0. In the case of / = 0, Eq. (fTSj) is reduced to K = 0, and the rest 
of equations can be integrated easily. In the case of / 7^ 0, it can be defined the variable 

^ <^ f = -. ,r, „r, , and Eqs. (fT3|) . (fH|l and (fT5]) can be rewritten as a 



y^fWTB ' V 1 - 

hyperbolic system of equations for F, A and b 



3. Conclusions 

Two basic results have been displayed: i) A theorem ensuring the existence of maximal slicings 
in any SSST. ii) A geometrical method to build up such slices by solving three decoupled first 
order quasi- linear partial differential equations p^ . (fH|) and p^ . The first result aims to fill 
a theoretical gap in the scientific literature. The second one tries to achieve an algorithmic 
procedure to obtain maximal slicings. An interesting by-product for Numerical Relativity of 
the approach presented in this paper has to do with the assessment of 3D codes written, as 
customary, in Cartesian coordinates. Let us consider two codes NCI and NC2 such that only 
NCI uses a gauge which is maximal. Hence, the evolution with code NC2 of any initial data 
admitting a spherically symmetric limit could be compared to the evolution produced by code 
NCI, by simply using our procedure to generate a SSST satisfying the maximal slicing condition. 
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